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Abstract 


Indirect estimates of gene flow from genetic data may substantially contribute to 
future analyses of population viability. However, depending on the model used, 
different assumptions have to be considered when interpreting the resulting 
estimates. 

We exemplify three indirect estimates of gene flow using a case study on allo- 
zyme divergence of the bush cricket Platycleis albopunctata in different land- 
scapes. The estimation of Nem from WRIGHT’s (1951) Fst and the graphical 
approach of SLATKIN (1981) produce comparable gene flow estimates. In con- 
trast, NEI's (1972) approach, which estimates the migration rate m from average 
genetic identities, produces values of Nam that are more than one order of magni- 
tude lower than the previous ones. 


Zusammenfassung 


Indirekte Methoden der Genflußschätzung, z.B. über die räumliche Verteilung 
alleler Varianten von Enzymgenorten, könnten in der Zukunft von grundlegender 
Bedeutung für die Abschätzung der Persistenz genetischer Variation in fragmen- 
tierten Landschaften sein. Sie basieren jedoch auf zahlreichen Annahmen, die bei 
der Interpretation der Ergebnisse berücksichtigt werden müssen. 

Am Beispiel einer populationsgenetischen Untersuchung an der Westlichen 
Beißschrecke Platycleis albopunctata werden drei indirekte Genfluß-Schätzver- 
fahren vorgestellt und ihre theoretischen Grundlagen erläutert. Die auf der F-Sta- 
tistik aufbauende Methode nach WRIGHT (1951) und die graphische Methode 
nach SLATKIN (1981) liefern Schätzungen der durchschnittlichen Anzahl immigrie- 
render Individuen pro Generation und Population (Nem), die in einer vergleich- 
baren Größenordnung liegen. Erstgenannte Methode überschätzt jedoch den 
Genfluß systematisch. Mit der Methode nach NEI! (1972) wurden wesentlich ge- 
ringere Genflußraten geschätzt, die auf einen unrealistisch hohen Isolationsgrad 
der Populationen in beiden Untersuchungsgebieten rückschließen lassen würde. 


Einleitung 


Das Arten- und Populationssterben in der Bundesrepublik Deutschland hält trotz 
intensiver Bemühungen zu ihrem Erhalt unvermindert an. Ein ständig steigender 
Flächenverbrauch und das stets dichter werdende Verkehrsnetz zerteilen die 
Landschaft in immer kleinere und isoliertere Lebensräume. Um dieser Entwick- 
lung Rechnung zu tragen, orientiert sich der moderne Artenschutz an der in einer 
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Landschaft ablaufenden Dynamik des Aussterbens und der Neugründung von 
Populationen (HENLE 1994). Über die Modellierung dieser Prozesse versucht man 
heute Wahrscheinlichkeitsaussagen zum Überleben von Populationen zu treffen 
(Populationsgefährdungsanalyse, PVA; siehe GILPIN & SOULE 1986, POETHKE et 
al. 1996). Von elementarer Bedeutung hierfür ist die realistische Abschätzung des 
Genflusses zwischen den Populationen. 


Eine Möglichkeit, Ausbreitungs- und Besiedlungsprozesse zu untersuchen, die 
dem lokalen Aussterben von Populationen entgegenwirken, ist die Analyse des 
Genflusses zwischen Populationen. Unter Genfluß versteht man das Einbringen 
von Genen einer Population in den Genpool einer anderen. Er setzt eine Ortsver- 
änderung von Individuen oder Gameten voraus. Der Begriff Genfluß wird daher 
häufig, aber fälschlicherweise, synonym mit dem Begriff Migration verwendet. 
Migration beschreibt nur die gerichtete Ortsveränderung von Individuen. Ihr Re- 
produktionserfolg bleibt hierbei unberücksichtigt. Die genetisch effektive Migra- 
tion, gemessen anhand des Genflusses, ist meist geringer als die tatsächlich 
stattfindende Bewegung zwischen den Populationen. Viele Individuen, die in eine 
andere Population immigrieren, haben dort aus verschiedenen Gründen keinen 
Fortpflanzungserfolg und treten deshalb genetisch nicht in Erscheinung. 

Im Falle von Tierpopulationen kann der Nachweis von Genfluß demnach nur ge- 
führt werden, wenn bewiesen wurde, daß (i) Individuen von einer Population in 
die andere gewandert und (ii) die eingewanderten Individuen reproduktiv erfolg- 
reich sind. Eine Bestimmung der Stärke des Genflusses ist nur möglich, wenn (iii) 
alle in eine Population immigrierenden Individuen erfaßt wurden und ihr repro- 
duktiver Erfolg bestimmt wurde. Mit populationsökologischen Methoden ist dies 
kaum zu bewerkstelligen. Folglich kommt der indirekten Bestimmung des Gen- 
flusses große Bedeutung zu. 


Aus der räumlichen Verteilung und Häufigkeit alleler Varianten einzelner Genorte 
lassen sich Genflußmaße schätzen. Im folgenden stellen wir verschiedene dieser 
Methoden vor und diskutieren anhand eines Fallbeispiels ihre Aussagekraft. Es 
soll dem Leser dadurch ermöglicht werden, auch ohne vertiefte Kenntnisse der 
Theorie, die Ergebnisse solcher Analysen - selbstverständlich mit aller Vorsicht - 
zu interpretieren. 


Indirekte Methoden der Genflußschätzung 


Es existieren mehrere Methoden der indirekten Genflußschätzung aus geneti- 

schen Daten: 

1. Schätzung der durchschnittlichen Anzahl der pro Generation und Population 
immigrierenden und sich erfolgreich reproduzierenden Individuen (Nəm, mit 
N= effektive Populationsgröße, und m = Migrationsrate) aus der genetischen 
Substrukturierung der Population (hierarchische F-Statistik), 

2. Schätzung von Nəm aufgrund des graphischen Vergleiches der durchschnitt- 
lichen Häufigkeit aller Allele, die in der gleichen Anzahl Populationen vor- 
kommen, mit unter bestimmten Migrationsraten simulierten Allelfrequenzen 
(SLATKIN 1981), 

3. Schätzung von m aus der genetischen Identität von Populationen (NEI 1972), 
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4. Schätzung von Nam mittels der durchschnittlichen Häufigkeit privater, d.h. nur 
in einer einzigen Subpopulation auftretender Allele (SLATKIN 1985) und 

5. Schätzung des Genflusses nach der Maximum Likelihood-Methode (BARTON et 
al. 1983). 

Die Anwendung der Methode der privaten Allele setzt die Kenntnis der Häufigkeit 

von mindestens 20 privaten Allelen voraus, um eine mit einem vernünftigen 

Fehlerbereich behaftete Schätzung zu erhalten (SLATKIN 1985). Die Maximum 

Likelihood-Methode liefert relativ ungenaue Schätzungen bei einer Anzahl an 

Genorten <40 (SLATKIN & BARTON 1989). Für die Anwendung beider Methoden ist 

daher schon die rein technische Voraussetzung nur äußerst selten gegeben, so 

daß wir uns im folgenden auf die ersten drei genannten Methoden beschränken. 


Genflußschätzung aus F;- 

Ein indirekter Weg der Genflußschätzung erfolgt über den von WRIGHT (1943) in 
die Populationsgenetik eingeführten Fixierungsindex Fs:. Populationen können in 
Subpopulationen aufgeteilt sein. Diese Unterteilung kann aufgrund eingeschränk- 
ten Genaustauschs zwischen den Populationen einen inzuchtähnlichen Effekt zur 
Folge haben, der sich anhand eines im Vergleich zu einer panmiktischen Ge- 
samtpopulation vorhandenen Defizites an heterozygoten Individuen nachweisen 
läßt. WRIGHTS (1943) F-Statistik mißt diesen Effekt der Abnahme des Hetero- 
zygotenanteils. Hierzu werden die drei Komplexitätsniveaus Individuum (1), Sub- 
population (S) und Totalpopulation (T) definiert. Fs , Fw und Fs, stellen nach 
WRIGHT (1951) die relativen Anteile an Heterozygoten einer Hierarchieebene an 
der erwarteten Heterozygotie der höheren Ebene dar. Für die Genflußschätzung 
ist nur der Fs: -Wert, auch Fixierungsindex genannt, interessant. Nur er wird 
daher hier ausführlicher beschrieben. 

Der Fixierungsindex Fs ist ein Maß für die Effekte der Populationsunterteilung in 
Subpopulationen. Er quantifiziert die Zunahme des Homozygotiegrades, bezogen 
auf die Totalpopulation, als Folge der Inzucht in den Subpopulationen. Die 
WRIGHTsche Formel für Fsr lautet: 


H, ist hierbei die erwartete Heterozygotie der Totalpopulation, H, ist die mittlere 


beobachtete Heterozygotie der Subpopulationen (vgl. auch die Darstellung in 
HARTL & CLARK 1989). 


Eine Abschätzung des Fsr -Wertes in Diallelsystemen erlaubt die Formel 


2 


> 


p-p) 
mit p = mittlere Frequenz eines der beiden Allele über alle Populationen und op? 
= Varianz der mittleren Allelfrequenz (WEIR & COCKERHAM 1984). 
Fs: wird genortweise berechnet und anschließend über alle Genorte gemittelt. 
Betragt Fs: = 0, befinden sich alle Subpopulationen an diesem Genort im HARDY- 
WEINBERG-Gleichgewicht und weisen die gleichen Allelfrequenzen auf, sie sind 
also panmiktisch. Bei Fs: = 1 sind alle Subpopulationen bezüglich eines Genortes 
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in einem Allel fixiert und in mindestens einer Subpopulation tritt ein von den übri- 
gen Populationen abweichendes Allel auf. 


Der Gebrauch von Fs: für die Schätzung des durchschnittlichen Genflusses 
basiert auf der Annahme, daß der Genfluß konstant ist und nur genetische Drift 
für den Unterschied zwischen Populationen verantwortlich ist. Als weitere, 
wesentliche Grundannahme geht man von selektionsneutralen Genorten aus; 
dies gilt auch für die anderen hier vorgestellten Schatzmethoden. Da jedoch Fsr 
als Mittelwert über die genortweise berechneten Fs,-Werte ermittelt wird, kann 
man anhand der zugrunde liegenden Einzelwerte Genorte mit stark abweichen- 
den Werten erkennen (Selektion, Interpretationsfehler, etc.) und diese bei der 
Mittelwertberechnung, wie von SLATKIN (1987) empfohlen, gegebenenfalls unbe- 
rücksichtigt lassen. 
Aus dem Fs,-Wert können die pro Generation migrierenden Individuen (N,m) ab- 
geschätzt werden. Im Inselmodell gilt nach WRIGHT (1951) für ein Drift-Genfluß- 
Gleichgewicht zwischen der Migrationsrate und der Verteilung genetischer Varia- 
tion zwischen den Populationen die Beziehung: 
| 
nme 

m+ | 

Durch Umformung ergibt sich als Genfluß: 


N,m= (+ - ) 
4\ Fa 


Wird Fs; anhand der Verteilung genetischer Varianten auf die Populationen er- 
mittelt, so läßt sich hieraus direkt Nem ableiten. WRIGHT (1978) setzt voraus, daß 
genetische Drift und Migration im Gleichgewicht stehen, d.h. die durch genetische 
Drift entstehenden Unterschiede werden durch Genfluß kompensiert. Nimmt der 
Wert für Nem einen viel kleineren Wert als 1 an, gelten die untersuchten Popu- 
lationen nach WRIGHT (1978) als voneinander isoliert. 

Im WRIGHTschen Inselmodell stellt genetische Drift die treibende Kraft bei der 
Differenzierung der Subpopulationen dar. Andere Kräfte, wie z.B. Mutation und 
Inzucht, treten in den Hintergrund (WRIGHT 1978, SLATKIN & VOELM 1991, WEIR & 
COCKERHAM 1984, 1993). Diese Annahmen gelten streng genommen jedoch nur 
für kleine Populationen. 

Neben der klassischen Methode nach WRIGHT (1943) gibt es weitere Algorithmen 
zur Berechnung von Fsr. WEIR & COCKERHAM (1984) vergleichen fünf solcher 
Methoden. Für Genorte mit mehr als zwei Allelen empfehlen sie ihr Maß @ für die 
Populationsunterteilung, ein aus den unterschiedlichen Komponenten genetischer 
Varianz (zwischen Populationen, zwischen Individuen einer Population, zwischen 
den Gameten eines Individuums) berechnetes Analogon zu WRIGHTS Fsr. 


Graphische Methode 

Als zweite Methode soll hier die graphische Methode nach SLATKIN (1981) vorge- 
stellt werden. Sie erlaubt ebenfalls die Abschätzung der durchschnittlichen An- 
zahl Migranten zwischen Populationen. Hierbei wird eine empirisch ermittelte 
Häufigkeitsverteilung der mittleren Häufigkeiten aller Allele, die in einem definier- 
ten Anteil an Populationen nachgewiesen wurden, mit unter vorgegebenen Gen- 
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flußraten ermittelten Modellrechnungen graphisch verglichen. Seine Computer- 
simulationen führte SLATKIN (1981) auf der Basis des "infinite allele models” von 
KIMURA & CROW (1964) durch. Es werden sowohl Mutationen berücksichtigt, die 
zu neuen Allelen führen, als auch solche, die den Verlust von Allelen nach sich 
ziehen. SLATKIN (1981) konnte mit seiner Simulation zeigen, daß die mittlere 
Allelfrequenz sowohl von der Selektionsstärke, als auch von der Mutationsrate 
relativ unabhängig ist. Sie hängt aber sehr stark vom Genfluß ab. Hieraus folger- 
te er, daß die mittlere Allelfrequenz für eine grobe Abschätzung des Genflusses 
in einer untergliederten Population verwendet werden kann. 

Bei der graphischen Darstellung werden die durchschnittlichen Allelfrequenzen 


aller Allele die in 'i' Populationen auftreten, gegen : aufgetragen. Dabei ist d die 


Anzahl aller untersuchten Populationen (vgl. Abb. 1). 


Schätzung der Migrationsrate m aus der genetischen Identität 
NEI (1972) stellte den Genfluß in Beziehung zur genetischen Identität /. Diese ist 
zwischen zwei Populationen an einem einzelnen Genort von NE! definiert als: 


2 xiy 


TS xP Dye | 

xj, yi sind die Häufigkeiten des i-ten Allels am j-ten Locus in den Populationen x 
und y. Da für eine sinnvolle Abschätzung der genetischen Differenzierung zwi- 
schen zwei Populationen mehrere Loci untersucht werden sollten, werden jeweils 
die arithmetischen Mittelwerte dieser Summen über alle untersuchten Loci ver- 
wendet. 

NEI (1972) nimmt an, daß die genetischen Identitäten durch Migration zwischen 
Populationen erhöht und durch Mutation erniedrigt werden. Die genetische Identi- 
tät und die zwischen zwei Populationen bestehende gegenseitigen Migrations- 
raten stehen in folgender Beziehung zueinander (NE! 1972): 

m, + m, 
m, +m, +2u 

Hierbei sind y die Mutationsrate und mı, m2 die gegenseitigen Migrationsraten 
zwischen zwei Populationen. Definiert man die mittlere Migrationsrate zwischen 
zwei Populationen als m- mem und setzt für u=10* pro Allozymgenort und Ge- 


neration ein (NEI 1972), so erhält man aus obiger Gleichung: 
1.2.10 ° 
1-I 

Für das Modell von Ne! müssen die Populationen groß sein, weil nur dann die 
Mutationsrate m den einzigen die genetische Divergenz beeinflussenden Faktor 
darstellt (Selektion wird vernachlässigt, da NEI von selektionsneutralen Allelen 
ausgeht). Natürliche Populationsgrößen sind immer kleiner als æ folglich tritt auch 
immer genetische Drift auf, womit die Modellannahmen in der Natur nie ganz er- 
füllt sind. Aus diesem Grund wird der Genfluß mit der NEischen Methode stets 
unterschätzt, da der Einfluß der Mutationsrate überbewertet wird. Für die Berech- 
nung von Nem muß die Migrationsrate lediglich mit einer unabhängig geschätzten 
mittleren Populationsgröße multipliziert werden. 
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Ergebnisse 


Im folgenden soll die Anwendung der voranstehend beschriebenen Methoden am 
Beispiel einer populationsgenetischen Untersuchung an der Westlichen Beif- 
schrecke, Platycleis albopunctata (GOEZE, 1778) dargestellt werden. Die Daten 
wurden an Populationen aus den Landschaftsräumen Haßberge (bei Zell am 
Main, Bayern), Mittelrheintal (zwischen Bingen und Neuwied, Rheinland-Pfalz) 
und Hammelburg (an der Fränkischen Saale, Bayern) gewonnen (SCHMELLER 





1995). 

- « - - Simulation (Nem = 0,25) 
— — Simulation (Nem = 2,5) 
— Zwischen Gebieten 
—— Haßberge 
—— Mittelrheintal 

i, 
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Abb. 1: Schätzung des Genflusses nach der graphischen Methode von SLATKIN 
(1981). N.= effektive Populationsgröße, m = Migrationsrate, p(i) = mittlere 
Häufigkeit aller Allele, die in i Populationen vorkommen, d = Gesamtzahl der 
untersuchten Populationen (vgl. auch Tab. 1). 


Fsr-Werte wurden hierbei nach WEIR & COCKERHAM (1984) für jeweils ein Gebiet 
berechnet. Für die Genflußschätzung zwischen den Gebieten wurde die hierar- 
chische F-Statistik nach WRIGHT (1978) angewendet. Für die Genflußschätzung 
nach NEI (1972) wurde m aus den genetischen Identitäten berechnet. Die Mittel- 
werte in Tab. 1 wurden aus den paarweisen Werten zwischen Populationen eines 
Gebietes (Haßberge, Mittelrheintal) bzw. zwischen Populationen unterschied- 
licher Gebiete berechnet. Zur Berechnung von Nem wurden im Rahmen popula- 
tionsökologischer Untersuchungen ermittelte Populationsgrößen N herangezogen 
(GOTTSCHALK, mündl. Mitt., sowie SCHMELLER 1995). Hierbei setzten wir 
näherungsweise N=N.. 
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Tab. 1: Ergebnisse der Schätzung von Nem mit verschiedenen Methoden. Mit der 
Methode nach Ne! (1972) kann nur die Migrationsrate m bestimmt werden. 
N.m wurde durch Multiplikation von m mit N unter der Annahme N, = N be- 
rechnet. In die Genflußschätzung zwischen Gebieten ging zusätzlich zu den 
Gebieten Haßberge und Mittelrheintal auch das Gebiet Hammelburg ein. 


| > Haßberge | Mittelrhein: 
Nm nach 
a7 | so | 20 
N | oro | 02 | 006 
sun | >25 | >25 |  =25 
weitere Parameter | | | 


Diskussion 

















Bevor Genflu&schatzungen im Kontext eines potentiell existierenden Popula- 
tionsverbundes diskutiert werden können, muß ihre Aussagekraft kritisch hinter- 
fragt werden. Eine indirekt über die räumliche Verteilung genetischer Varianten 
geschätzte Migrationsrate Nəm berücksichtigt auch Geschehnisse, die zeitlich 
weit zurückliegen können. Diese Form der Genflußbestimmung stellt somit „ein 
Integral über Raum und Zeit” dar (DUNLEY & CROFT 1994, VEITH & SEITZ 1995). 
Bei früheren Massenwanderungen können sich z.B. die Individuen unterschied- 
licher Populationen vermischt haben, während sie anschließend für mehrere 
Generationen nicht mehr migrierten. Der geschätzte Genfluß wird dann über 
sämtliche Generationen gemittelt, d.h. man errechnet auch für solche Generatio- 
nen einen Genfluß, in denen de facto keiner stattgefunden hat. Wie gut dieses 
„Gedächtnis“ funktioniert, hängt nicht zuletzt von der Stärke der genetischen Drift, 
d.h. von den lokalen Populationsgrößen ab. 

Ferner kann ein doppelter Zusammenhang zwischen der Populationsgröße und 
der Anzahl migrierender Individuen bestehen. Einerseits ist bei großen Popula- 
tionen auch bei konstanten Migrationsraten das Produkt Nem größer, andererseits 
kann die Migrationswahrscheinlichkeit gegenüber kleinen Populationen zusätzlich 
erhöht sein (dichteregulierte Migrationswahrscheinlichkeit; POETHKE, münd. Mitt.). 
Für die von uns untersuchten Populationen liegen leider keine langjährigen Erhe- 
bungen zu Populationsgrößen vor. Es ist aber bekannt, daß Populationen von 
P. albopunctata sehr stark fluktuieren können (bis Faktor 10 in aufeinander- 
folgenden Jahren; GOTTSCHALK mündl. Mitt.). 

Die von uns über die Migrationsrate m (NEische Methode) geschätzten Migran- 
tenzahlen weichen stark von den beiden anderen Werten ab (Tab. 1). Das 
NEische Migrationsmodell zeichnet sich durch einen wichtigen Unterschied zu 
den Grundannahmen des Inselmodells aus: NEI sieht nicht die genetische Drift 
als treibende Kraft der Differenzierung, sondern die Mutation. Dieses Modell 
sollte daher besser geeignet sein für Genflußschätzungen zwischen großen Po- 
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pulationen. Zudem geht die Mutationsrate y in die Schätzung ein, über deren 
Höhe bei verschiedenen Arten und Umweltbedingungen noch wenig bekannt ist. 
Der von uns eingesetzte Wert wurde empirisch bei Fruchtfliegen (Drosophila) 
ermittelt und bezieht sich auf die Detektion neuer, elektrophoretisch nachweis- 
barer Allozymvarianten (z.B. FUTUYMA 1990, HARTL & CLARK 1989). 

Unter der Annahme, daß sowohl die Mutationsrate als auch die eingesetzte effek- 
tive Populationsgröße realistisch sind, ergibt sich aus dem von uns geschätzten 
Wert von Nem = 0,10 in den Haßbergen ein Fs: von ca. 0,7. Ein solch hoher Wert 
würde jedoch für hochgradige Isolation zwischen den Subpopulationen sprechen, 
eine angesichts der gefundenen Allelfrequenzen unrealistische Annahme. Zieht 
man zudem in Betracht, daß die Annahme N = N, ebenfalls unrealistisch ist, da 
insbesondere bei stark fluktuierenden Insektenpopulationen N.<<N ist (bei 
schnell wachsenden Populationen wird N, maßgeblich durch die zwischenzeitlich 
kleinste effektive Populationsgröße bestimmt; siehe hierzu auch HARTL & CLARK 
1989), läge ein aufgrund einer realistischeren effektiven Populationsgröße 
berechneter Wert für Nem noch niedriger, womit zwangsläufig Fs, noch näher an 
den Wert 1 heranrücken würde. 


Zieht man die Genflußmodelle von NEI und WRIGHT in seine Überlegungen ein, 
so wird man feststellen, daß sie Grenzfälle darstellen. Das WRIGHTsche Modell 
überschätzt den Genfluß zwischen großen Populationen, weil die genetische Drift 
in diesen nicht mehr die dominierende Kraft für die genetische Diversität ist. Zu- 
dem führt die Verwendung von WEIR & COCKERHAMS (1984) 8 ebenfalls zu einer 
systematischen Überschätzung von Nam. Das NEische Modell unterschätzt den 
Genfluß, wenn die Populationen nicht unendlich groß sind, da der Einfluß der 
Mutationsrate überschätzt wird. In der Natur wird keiner dieser Grenzfälle er- 
reicht, sondern es wirken immer sowohl die Mutation, wie auch die genetische 
Drift, je nach Populationsgröße mit unterschiedlicher Bedeutung. Simulations- 
experimente (SLATKIN & BARTON 1989, SLATKIN 1993) sowie unsere Überlegun- 
gen im voranstehenden Abschnitt legen jedoch nahe, daß eine Genflußschätzung 
über Fs: die realistischeren Werte liefert. Die Abschätzung von Nam mittels der 
graphischen Methode führt zu Ergebnissen, die den mit Hilfe der F-Statistik 
berechneten Werten vergleichbar sind. 


Der Wert Nem, als „migrierende Individuen pro Generation” bezeichnet, ist kein 
exaktes Maß für die wanderenden Tiere, sondern ein Relativwert, der sich aus 
historischen und gegenwärtigen Ereignissen ergibt. N.m sollte in bezug auf 
gegenwärtige Migrationsbewegungen nicht überbewertet werden. Es empfiehlt 
sich daher, die Historie einer Art im untersuchten Landschaftsraum zu berück- 
sichtigen, um zu geeigneten Schlußfolgerungen zu gelangen. Innerhalb eines 
Gebietes bieten Genflußschätzungen allerdings eine gute Vergleichsmöglichkeit 
zwischen Arten (vgl. VEITH et al. 1996). 
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